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Abstract 

Starting from the conventional electron-hole Hamiltonian 7i e h, we derive an 
effective Hamiltonian 7i\ s for Is excitons with spin degrees of freedom. The 
Hamiltonian describes optical processes close to the exciton resonance for 
the case of weak excitation. We show that straightforward bosonization of 
TL e h does not give the correct form of H\ s , which we obtain by a projection 
onto the subspace spanned by the Is excitons. The resulting relaxation and 
renormalization terms generate an interaction between excitons with opposite 
spin. Moreover, exciton-exciton repulsive interaction is greatly reduced by the 
renormalization. The agreement of the present theory with the experiment 
supports the validity of the description of a fermionic system by bosonic fields 
in two dimensions. 

PACS numbers: 71.35.-y, 71.10.-w, 78.66.-w 
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Properties of a fermionic system can sometimes be described by bosonic fields, by which 
theoretical analyses may be greatly simplified. The most successful example is the bosoniza- 
tion of one- dimensional (d — 1) conductors JT[, where success is due to the specific feature 
of the pair spectrum |2[], i.e., the low-energy pair spectra consist of discrete branches in 
the energy versus momentum plane. For d > 2 conductors, discrete branches overlap with 
continuous spectra, hence the bosonization is nontrivial and still in progress ||. 

For insulating solids, on the other hand, discrete branches of excitons are separated 
from the continuum for any d. From this point of view, it has been suggested that a 
useful bosonic theory may be constructed if one focuses on exciton states, even for d > 2 
J§^7|]. However, the validity of the bosonic description of excitons is nontrivial, because the 
binding energy (in, e.g., GaAs) is comparable to other relevant energies ||. Under an optical 
excitation, excitons and free electron- hole (e-h) pairs (continuous spectra) will be created. 
As the excitation intensity (and thus the e-h density) is increased, the fermionic nature of 
the system becomes more important, and bosonization requires more bosonic fields. Such 
a strong-excitation regime has been successfully analyzed without the use of bosonization 
|]| |TB[. In a weak-excitation regime, on the other hand, it is expected that the system is 



well described by a small number of bosonic fields. If this is the case, the bosonic theory 
will provide a powerful theoretical tool as well as a transparent physical view, as in the case 
of the d — 1 conductors [p]] . 

To demonstrate the effectiveness of the bosonized theory, optical experiments may be 
more convenient than the electron transport experiments, because one can easily produce 
and detect two or more light beams, and obtain rich information from responses to the 
multi beams. Moreover, one can easily control the polarization of each individual light 
beam, which gives more detailed information. Recently, by controlling the polarizations of 
two light beams, Kuwata-Gonokami et al. demonstrated experimentally that the polariton- 
polariton scatterings in a quantum well (QW) in an optical micro cavity are well described 
by a phenomenological bosonic Hamiltonian in the weak-excitation regime pJTl]. Their 



experiment strongly indicated the validity of a bosonic description for d = 2. However, no 



theoretical studies were reported which derive their phenomenological Hamiltonian from a 
microscopic fermionic theory. 

In this letter, we derive an effective Hamiltonian of excitons from the conventional e-h 
Hamiltonian. In particular, we will show that the derivation of an effective bosonic theory 
is nontrivial and not straightforward because a direct bosonization does not give correct 
results. Our calculation and the agreement with the experimental data f| directly prove 
the validity and relevance of a bosonic description for d = 2. 

Model — We start from a Hamiltonian 7i e h for an interacting electron-hole system in a 
QW. 

Here, V(r) denotes the Coulomb potential, which behaves in a QW of width L as V(r) ~ 
e 2 /er for |r| > L, where e is the static dielectric constant, and V(r) ~ constant for |r| < L. 
We simplify the calculation by taking the limit L — > wherever the singularity at r = 
is irrelevant. In eq. (JJ), ^ e ^)(x) is the field operator of an electron (hole), {i,i'} = {e, h}, 
z e {h) = x = (r e (h), Je(h))> and Jdx = Ej*/d 2 r». The index J| (h) denotes the z- 

component of the total angular momentum (which is a good quantum number), where we 
take the z-axis to be in the direction normal to the QW layers. The is defined as 
— 1 times J z of a valence band electron. In a GaAs QW, e.g., = ±3/2 for a heavy 
hole, and J| = ±1/2 p. A photon with J^ h = ±1(— 1) creates an electron-hole pair with 
J z e = -1/2 (±1/2) for the electron and with J% = +3/2 (-3/2) for the hole. The exciton 
states are labelled by indices q, v and S, where q is the momentum of the center-of-mass 
motion, v denotes the set of quantum numbers for the relative motion {y = Is, 2p + , 2p_, ■ ■ •) 
||, and S denotes combinations of J| and J%. Since S is related to the total angular 
momentum, we call S the "spin" index. 

Since all states excited by photons are electrically neutral, we can confine ourselves to the 
charge-neutral sector. We first transform 7i e h to a Hamiltonian H. of interacting excitons. 
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We then make a projection onto the subspace spanned by Is excitons to obtain an effective 
theory which is described by a renormalized Hamiltonian 7Y ls of Is excitons and a relaxation 
f. 

Transformation of 7i e h into Ti. — We first rewrite 7i e h as a function of exciton operators 

r _ f A 2 12 1 (■ m e r e + m h r h \ 
\„s = E yd r e d r™ exp • j 



7Z /Z 

x^(r e - r h )<S|J*, J^)4(r e , J*)4(r fc , (2) 

Here, <p v {v) is a wavefunction for the e-h relative motion, (S\J*,Jf l ) the Clebsch-Gordan 
(CG) coefficient, Q the QW area, and M = m e + nih. In the following, we assume < m e <C 
m/j. To define S uniquely, we choose S = +,—,a, and /3, for \ J%, J|) = | — 1/2, +3/2), 
1 + 1/2, -3/2), | + l/2, +3/2), and [ — 1/2, -3/2), respectively, in a GaAs QW @. An exciton 
with S — ± couples to circularly polarized light with = ±1, whereas excitons with S = a 
or (3 are dipole inactive. Then, assuming that the exciton density is low, we regard & qt ,s' s 



as genuine boson operators |15j and rewrite H e h in the charge-neutral sector as H e h — * TC = 
HP + H, int , where 7i° and 7i mt are the free and interaction parts of excitons, respectively. 
Since the transferred momentum in the exciton scattering processes is fairly small, and of the 
order of the photon momentum, the direct and the double fermionic exchange interactions 
are negligible and the momentum dependence of the exchange interaction can be omitted. 
Hence, the general form of Ti mt is given by 

n mt = E S E (s^J^is^'jnis^'JOis^j*') 

kk'q{S}M ZU JiJWJZ 

X ^k+qvi Si kk'-q^kk'^sA^Ss • (3) 

Using this formula, we express Ti as 

H = Ti° + H ls + H' ls + Hothers, (4) 

where Tif s and 7i' ls include the v — Is operators only, and Ti ot hers denotes the remaining 
terms. The 7^f s consists only of the operators with S = ±, whereas W ls consists of terms of 
5 = a and /? operators, including cross terms with S = ± operators: 



H ts= ^EE 6 k+qs4'-qsVs 6 kS» ( 5 ) 
Zi L S=± kk'q 



n ' ls = Ti ^ 

" kk'q 



5E ( 2^k+q5^k'-q5^k'5^k5 + ^L+q +^k'-q 5^k' 5^k + 



.5=a,/3 

fo k+q-^-q5 & k' A-) + ( 6 k+q+&k'-q- 6 k'A/3 + H - C - 
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where we have written & q i s s = ^ q s, an d the effective interaction strength U = U v= \ s is 
expressed as 

u = 2v^E v(p - p') [|^ s (p)|VL(p)^(p') - I^(p)I 2 IvMp')I 2 ] • (7) 

p,p' 

Here, V"(p) and <p u (p) are the Fourier transforms of V(r) and ip u (r), respectively. For 
d — 2, and in the limit of L — > 0, eq. (|^) is evaluated as Z7 = 1.52aQ^ |I7| , where 
a = eh 2 M/e 2 m e mh and E 1 ^ = 2e 2 /eao- 

It is tempting to take 7if s + Ti.' ls as an effective interaction Hamiltonian Ti^J for Is exci- 
tons, by simply dropping 7i thers- However, as we will show below, this is not an appropriate 
approximation of Tiig , because important terms which describe the interactions between 



excitons with S = + and with S — — are absent in Hf s + 7i' ls |18 |. 

Projection — The evolution of the density operator p in the charge neutral sector obeys 
the von Neumann equation, dp/dt = (l/ih)[H, p\. To derive the correct His , we make a 
projection onto the subspace spanned by the Is excitons. We then obtain the equation for the 
reduced density operator p ls , as dp\ s /dt = (l/ifi)[Hi 8 , Pis] + ^Pis- Here, Hi s describes the 
unitary evolution of p~i s , and T the relaxation operator. We can calculate optical responses 
using His, neglecting T, when (a) the photon energy is close to the energy of the Is excitons, 
(b) the photoexcitation is weak, and (c) relaxation processes do not play an important role 
during the photon scattering processes. The reason for (a) is that 2p and higher excitons 
have been projected out. The condition (b) is due to the fact that the deviation from the 
Bose statistics of operators b^s and b^ vS becomes non-negligible when the excitation is 
strong ||19|| . In the experiment described in ref. ||, the three conditions have been satisfied. 
To satisfy condition (c), both the resonance of Is excitons and the strong coupling to the 



radiation field in a high Q micro cavity are crucial p0[ . Otherwise, we must take account 
of T in the calculations of optical responses. 
The renormalization procedure yields pT | 



(8) 



where 7Y° S is the Hamiltonian of free Is excitons, and Hf s and 7i' ls include the v = Is 
operators only. The T-Cf s consists only of operators with 5 = ±, whereas H' ls consists of 
terms of 5 = a and (3 operators, including cross terms with 5 = ± operators. Since TC' ls 
includes dipole inactive Is excitons (5 = a, (3), it does not contribute to the optical response 
in its lowest order [14,22]. The most important terms are therefore included in Ti.f s , which 
is evaluated, to the second order in the exciton-exciton interactions, as 



± U-U'^i t U' 



2tl 
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S=± kk'q 

where U' is a positive constant which arises from the renormalization of higher exciton states 
(u = 2p + ,2p_, •••): 



K ^ ls 2 (£„ + K 2 /2M) - 2£ ls 



Z ^(p - p + K ) h^L(p)^L(p')^(p)^(p') 



p,p' 



+ 2#,(p)#,(p - K)^(p)^(p') -y$L(p)v5L(p')v5,(p - K)^(p' + K) 



(10) 



Comparing the right-hand side of eq. @ with that of eq. (|5|), we find that the coefficient of 
the first term is renormalized as U — > U — U', and that a second term is generated which 
leads to an interaction between the 5 = 4- and — excitons. That is, the renormalization 
of higher exciton states results in the renormalized Hamiltonian Tifg, which differs, both 
quantitatively and qualitatively, from the bare Hamiltonian 7if s . 

We argue that the correct form of the effective Hamiltonian for Is excitons is the renor- 
malized one, i.e., TixJ = Hf s + H' ls . In fact, the interaction between the 5 = 4- and — 
excitons in 7t£ 1s , which is absent in 7i\ s , has been clearly observed experimentally in refs. 
H and P3| . Kuwata-Gonokami et al. expressed this interaction as an interaction term 
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(whose coupling constant is W) in the phenomenological Hamiltonian, which also has the 
interaction term (whose coupling constant is R) of excitons with parallel spins. The phe- 
nomenological Hamiltonian has the same form as 7Y ls , the dipole active part of TCu- This 
is quite reasonable because the other part Ti[ s , which is dipole inactive, should be invisible 
in low-order optical experiments P2 |. We can, therefore, identify the parameters R and W 



of the phenomenological Hamiltonian || as 

R=(1.52a 2 E b x -U')/(2Q), (11) 

w = -u'/n. (12) 

The value of U', as given by eq. ([Top, depends on the material parameters such as M and 
e, and hence it is different for different materials. It also depends on the QW parameter L. 
Moreover, when imperfections in the QW are non- negligible, the formula for U' should be 
modified accordingly. Therefore, even for the same material, the values of R and W could 
vary slightly from sample to sample, which seems to be consistent with recent experimental 
results [J2tJ. Note, however, that the existence of both terms of Hf s is independent of such 
details. 

We here estimate the typical value of U' as follows. The K-summation in eq. ( ]TD| ) is 
cutoff for K > Cl/L (through the cutoff of V) and/or for K > C ao /a (through (p v \ where 
Cl and C ao are cutoff parameters of the order of unity. For the case of the QW sample of 
ref. 0, L m ciq, hence we may cutoff the K-summation for K > C/clq, where C is of the order 
of unity. For the v summation, we may consider v = 2p± states only, because higher exciton 
states give much smaller overlap integrals. These approximations yield U' ~ 16.5C 2 alE^.. 
On the other hand, ref. |J reported the ratio R : W as 1 : —15. From eqs. (|TTD and (|12~D , we 
find that this ratio is reproduced by the present theory when the cutoff parameter C ~ 0.3, 
which is consistent with the requirement that C is of the order of unity. Considering that 
the values of R and W vary slightly from sample to sample [p0[] , the agreement seems 
satisfactory. Note that such a small value of R is due to the renormalization of U — > U — U'. 
Once the agreement of Ti.f s with the phenomenological Hamiltonian is thus established, the 



agreement with the experiment follows. That is, lowest-order perturbational calculations for 
the polariton-polariton scattering amplitudes agree with the experiment PJI^|. 

Discussions and remarks — It was conjectured |J that a "biexciton effect" would be 
the origin of the " W term" , the interaction between S — + and — excitons. However, this 
argument is misleading. The biexciton state is formed essentially from the mixing of two 
Is states having different centers. For examples, in the case of a hydrogen molecule, the 
mixing yields the bonding and antibonding states, (l/\/2) (cLcL ± c^^-^hi^h^,^)) . Here, 
c^ 2 ) creates an electron in the Is state located at nucleus 1(2), and h\, 2 \ creates the nucleus. 
In the case of excitons with J| = ±1/2 and = ±3/2 the corresponding states are 
(l/v^^+fr- ± &Q&g]|0), where we have not shown the k-dependence in order to focus on 
the S'-dependence. The bonding state (— sign for a positive coupling constant) has a lower 
energy and is called a biexciton. This energy splitting between the bonding and antibonding 
states is induced by the interaction of the form of b\_bLb a bp + H.c, which is included in 7i' ls 
(or, before the renormalization, in H/ ls of eq. @) On the other hand, the W term 



lowers the energies of both states by the same amount, hence does not play a central role in 
the formation of the biexciton state. The most important effect of the W term is to lower 
the energy of fe+&!_|0), relative to those of 6+6+|0) and 6L&L|0), and this effect was detected 
experimentally [f|. In short, in the framework of the present theory, 7i' ls lowers the energy 
of the bonding (biexciton) state relative to that of the antibonding state, and thus is crucial 
for the formation of the biexciton state, whereas the W term lowers the energy of both 
bonding and antibonding states. 

Note that 7ii s is not positive definite to the fourth order in the exciton operators. The 
stability of the system should be preserved by higher order terms. In general situations, 
properties of a system described by such a Hamiltonian should not be analyzed by a per- 
turbation theory based on the vacuum of the free part. Nevertheless, we can use such 
a perturbation theory in our case, because our exciton theory has the built-in constraint 
that the ground state is the state with no excitons, i.e., the vacuum of 7l^ s . The effective 
Hamiltonian 7ii s together with this constraint constitutes a consistent theory, which justifies 
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the low-order perturbation theory based on the given vacuum, if the optical excitation is 
sufficiently weak. 

We have used a low-order perturbation theory to derive 7ii s . However, this does not 
imply a total neglect of higher order terms, because we have calculated a Hamiltonian 
rather than observables. In fact, a systematic summing up of higher order terms is already 
incorporated in our theory if one calculates higher order scattering amplitudes, e.g., by 
writing the Bethe-Salpeter equation and using 7ii s . 

Finally, we discuss the relation between the fermionic theories PHl3| and our bosonic 



theory. The Hartree-Fock (HF) factorization treatment of the semiconductor Bloch equa- 
tions |l(J can not produce the interaction between the S = + and — excitons. The HF 



theory, therefore, corresponds to 7if s , eq. flj^). It was argued in refs. [p7T| P~3T that the in- 
teractions of an exciton with higher states (including free carriers) are important, and that 
the interactions result in the energy shift, the excitation-induced dephasing (EID), and the 
"biexcitonic correlations". In the bosonic theory in the form of eq. (^), these effects are 
included in H' ls and 7i thers- After the projection is made, the relation is roughly as follows. 
The renormalized Hamiltonian Ti.f s , eq. (0), would include the HF term and a part of the 
"biexcitonic correlation." The EID may be described by both F and 7~C[ S . Another part of 
the "biexcitonic correlation" would also be included in 7i' ls . The present theory thus helps 
to bridge the gap between the bosonic theories 0-0] and the fermionic theories [p|-[T3|] of e-h 
systems. However, more detailed comparisons will be a subject of future studies. 

Summary and Conclusions — In this letter, starting from the conventional electron-hole 
Hamiltonian (eq. (1)), we have derived the effective Hamiltonian 7ii s for the Is excitons 
with the spin degrees of freedom using the bosonic exciton operators in two dimensions (eq. 
(2)). It is found that the renormalization associated with the projection onto the Is exciton 
space is crucial, which leads to the generation of the attractive interaction between excitons 
with opposite spins (eq. (12)), and to the large reduction of the repulsive interaction between 
excitons with parallel spins (eq. (11)). Such a drastic modification of the interactions was 
absent in the previous theory without the renormalization procedure. The present theory 



is valid for systems that satisfy the following conditions: (i) excitation is weak, (ii) the Is 
excitons play a crucial role, and (iii) the exciton relaxation process is less important due to, 
for example, the micro cavity of a high Q-value. This effective Hamiltonian provides the 
microscopic foundation of the phenomenology proposed in ref. ||. The agreement of the 
present theory with the experiment supports the validity of the description of a fermionic 
system by bosonic fields in two dimensions. 

Helpful discussions with Professor Kuwata-Gonokami and Dr. Suzuura are acknowledged. 



10 



REFERENCES 



* Electronic address: inoue@ASone.cu-tokyo.ac.jp 

** Present address: 1. Inst. Theor. Physik, Jungiusstr. 9, D-20355 Hamburg, Germany 
[1] M. Stone: Bosonization (World Scientific, Singapore, 1994). 
[2] F. D. M. Haldane: J. Phys. C14 (1981) 170. 

[3] P. Kopietz: Bosonization of Interacting Fermions in Arbitrary Dimensions (Springer, 
Berlin, 1997). 

[4] E. Hanamura and H. Haug: Phys. Rep. 33 (1977) 209. 

[5] M. Kuwata-Gonokami, S. Inouye, H. Suzuura, M. Shirane, R. Shimano, T. Someya and 
H. Sakaki, Phys. Rev. Lett. 79, 1341 (1997). 

[6] V. M. Axt and S. Mukamel: Rev. Mod. Phys. 70 (1998) 145. 

[7] A. L. Ivanov, H. Haug and L. V. Keldysh: Phys. Rep. 296 (1998) 237 and references 
therein. 

[8] This should be contrasted with liquid Helium 4: there is no doubt that an Helium 
4 atom can be described by a boson operator when the superfluidity is studied, be- 
cause its nuclear binding energy is much larger than typical energy scales of superfluid 
experiments. 

[9] See, for example, H. Haug and S. W. Koch: Quantum Theory of the Optical and Electric 
Properties of Semiconductors (World Scientific, Singapore, 1994) 3rd ed. 

[10] M. Lindberg and S. W. Koch: Phys. Rev. B38 (1988) 3342. 

[11] Y. Z. Hu, R. Binder, S. W. Koch, S. T. Cundiff, H. Wang and D. G. Steel, Phys. Rev. 
B49, 14382 (1994). 

[12] T. Rappen, U. B. Peter, M. Wegener and W. Schafer: Phys. Rev. B49 (1994) 10774. 

11 



[13] W. Schafer, D. S. Kim, J. Shah, T. C. Damen, J. E. Cunningham, K. W. Goossen, L. 
N. Pfeiffer and K. Kdhler, Phys. Rev. B53, 16429 (1996). 

[14] It should be stressed that fairly good agreement is obtained in ref. between the 
experimental results and the lowest-order calculations based on the phenomenological 
Hamiltonian. 

[15] H. Haug and S. Schmitt-Rink: Prog. Quant. Electr. 9 (1984) 3. 



[16] We have obtained a similar result for the case of Jh = ±1/2 pi . 



[17] S. Schmitt-Rink, D. S. Chemla and D. A. B. Miller: Adv. Phys.39 (1989) 89. 

[18] Ivanov and Haug obtained this interaction without performing the projection (Phys. 
Rev. B48 (1993) 1490). However, they dropped CG coefficients of eq. (|3|), and deter- 
mined the spin dependence using another theory of a four-body problem. Hence their 
derivation is not fully microscopic. 

[19] M. Hawton and D. Nelson: Phys. Rev. B57 (1998) 4000. 

[20] M. Kuwata-Gonokami: private communication. 

[21] J. Inoue, T. Brandes and A. Shimizu (unpublished). 

[22] However, TC' ls , of course, becomes important when one measures optical processes in 
which higher-order terms are involved. 

[23] K. Bott, O. Heller, D. Bennhardt, S. T. Cundiff, P. Thomas, E. J. Mayer, G. O. Smith, 
R. Eccleston, J. Kuhl, and K. Ploog, Phys. Rev. B48, 17418 (1993). 

[24] A typical bound state in an interacting boson model is the eigenstate \(f>) = (a^tf — 
sign(g)cW)\0) oiH = v(da + tfb + Jc + tfd)+g(aWcd + (tfctba). The eigenenergy is 
2uj — \g\. Note that the existence of the bound state is independent of the sign of the 
interaction g. 



12 



